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Suppose that a d-dimensional domain is filled with a gas of (in general, interacting) diffusive par¬ 
ticles with density no- A particle is absorbed whenever it reaches the domain boundary. Employing 
macroscopic fluctuation theory, we evaluate the probability V that no particles are absorbed during 
a long time T. We argue that the most likely gas density profile, conditional on this event, is sta¬ 
tionary throughout most of the time T. As a result, V decays exponentially with T for a whole class 
of interacting diffusive gases in any dimension. For d = 1 the stationary gas density profile and V 
can be found analytically. In higher dimensions we focus on the simple symmetric exclusion process 
(SSEP) and show that — ln"P ~ DoTL d ~ 2 s(no), where Do is the gas diffusivity, and L is the linear 
size of the system. We calculate the rescaled action s(no) for d = 1, for rectangular domains in 
d = 2, and for spherical domains. Near close packing of the SSEP s(no) can be found analytically 
for domains of any shape and in any dimension. 

PACS numbers: 05.40.-a, 02.50.-r 


I. INTRODUCTION 

Diffusive lattice gases serve as useful simplified mod¬ 
els of many stochastic spatio-temporal systems in nature. 
Among them are diffusion-controlled chemical reactions: 
reactions which occur quickly once the diffusing reagent 
particles “find” each other in space. A simple but amaz¬ 
ingly rich model of this process, due to Smoluchowski 
[I], treats one of the two reacting species as an immo¬ 
bile large-size minority. The Smoluchowsi’s model allows 
one to calculate the expected reaction rate (that is, the 
expected rate of absorption of a random walker by a tar¬ 
get). Statistics of fluctuations of this rate have been the 
subject of numerous studies When the majority 

molecules are treated as noninteracting random walkers, 
the calculation of the effective reaction rate and its fluctu¬ 
ation statistics boils down to calculating a single-particle 
probability. Recently some progress has been also made 
in the situation when the diffusing particles interact with 
each other [7| d] . 

Here we extend this line of work by putting the walkers 
inside a domain which boundary is the “target”. Suppose 
that a gas of diffusing and, in general, interacting parti¬ 
cles with density n o fills a d-dimensional domain Cl. Each 
particle is absorbed immediately whenever it reaches the 
domain boundary dCl. This simple setting is a carica¬ 
ture of a host of processes inside the living cell, where a 
molecule needs to reach the cell membrane )9f. We as¬ 
sume that, on macroscopic length and time scales, the 
average gas density inside the domain, n(x, t), evolves 
according to a diffusion equation, 

d t n = V ■ [D(n)Vn] (1) 

with diffusivity D{n) that may depend on n. The bound¬ 
ary condition at the absorbing domain boundary is 

n(x € dCl, t) = 0. (2) 

Solving Eqs. © and © for a given initial condition, such 
as n(x, t = 0) = no, and calculating the diffusion flux into 


the domain boundary, one can find the expected number 
N(T) of absorbed particles during time T. In individual 
realizations of the underlying microscopic stochastic pro¬ 
cess, the number of absorbed particles fluctuates around 
N(T), and it is interesting to determine the fluctuation 
statistics. In this work we will deal with an extreme limit 
of this statistics, corresponding to the survival probabil¬ 
ity V(T): the probability that not a single particle hit 
the domain boundary by time T which is long compared 
to the characteristic diffusion time through the domain. 
For non-interacting diffusing particles (we will call them 
Random Walkers, or RWs), D{n) = D 0 = const. In this 
case one obtains @ 

— IuPrw — uqDqTV p\, (3) 


where V is the domain volume, and pi is the lowest pos¬ 
itive eigenvalue of the eigenvalue problem \7 2 u + p 2 u = 0 
for the Laplace’s operator inside the domain with the 
boundary condition u(x £ dCl) = 0. For a d-dimensional 
sphere of radius R , Eq. © yields the well-known results 

-\nV RW ~n 0 D 0 TR d - 2 f d , (4) 

where 

d= 1, (5) 

d= 2, (6) 

d = 3, (7) 

and Z\ = 2.4048 ... is the first positive root of the Bessel 

function Jq{z). 

The exponential decay of 'Prw with time T, as de¬ 
scribed by Eqs. © and ©, reflects the fact that, at long 
times, the single-particle survival probability decays ex¬ 
ponentially with time. Indeed, for a single RW, the sur¬ 
vival probability distribution p\ (x, t) inside the domain 
obeys the diffusion equation d t p\ = D 0 X7 2 pi with the 
absorbing boundary condition pi(x £ dCl,t) = 0 and 
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a delta-function initial condition Q. Finding pi(x, t) 
and integrating it over the domain, one obtains the 
single-particle survival probability as a function of time. 
Its long-time asymptotic describes an exponential decay 
with the decay rate corresponding to the lowest positive 
eigenvalue ^1 of the Laplace’s operator. The gas survival 
probability Vew is given by the product of the survival 
probabilities of all independent particles inside the do¬ 
main. What is left to arrive at Eq. © is to go over to 
the continuum limit by replacing a discrete sum in the 
exponent of Prw by an integral. For completeness, this 
procedure is presented in Appendix [A] As a result, the 
exponential long-time decay of the survival probability 
holds, for the independent RWs, in all spatial dimensions, 
as evidenced by Eq. ©• 

An important finding that we report here is that the 
exponential-in-time decay of the survival probability V 
holds when the diffusing particles interact with each 
other, and the single-particle picture breaks down. This 
non-trivial result is a consequence of the fact that the 
optimal gas density history, conditional on the long-time 
survival of all particles, is almost stationary, in any spa¬ 
tial dimension. We show it here by employing the Macro¬ 
scopic Fluctuation Theory (MFT) [10]. For d = 1 the 
stationary MFT problem is soluble for a whole family of 
interacting gases. In higher dimensions the solution is in 
general unavailable. Here we focus on a gas of particles 
interacting via exclusion, so as to describe e.g. diffusion- 
controlled chemical reactions in a crowded environment 
of a living cell Q Specifically, we will study the sym¬ 
metric simple exclusion process (SSEP). In the SSEP 
each particle can hop to a neighboring lattice site if that 
site is vacant. If it is occupied, the move is not allowed 
[hi . For the SSEP we obtain — InP ~ D^TL d ~ 2 s(?ro), 
where L is the linear size of the domain [13|. We calcu¬ 
late the rescaled action s(no) for several domain shapes 
and in different dimensions. As we show, s(n o) increases 
with the density no faster than linearly, see Fig. [5] for 
d = 1, Fig. 0 for a rectangle in d = 2, and Fig. [9] for a 
sphere in d = 3. In the MFT formalism s(no) diverges as 
no approaches unity, but this divergence is cured when 
no is sufficiently close to unity where the MFT breaks 
down. 

The interior survival problem, considered here, has an 
exterior analog that is known by the name “survival of 
the target”. In the exterior problem a gas of particles 
surrounds an absorbing domain from outside. As in the 
interior problem, one is interested in the probability that 
no particle hits the domain during time T. For the RWs 
the exterior problem was extensively studied in the past 
[!H23. Recently, the theory has been extended to in¬ 
teracting diffusive gases: for the survival probability Q 
and in the more general context of full absorption statis¬ 
tics j[8j. As our present work, Refs. 00 employed the 
MFT formalism. In contrast to the interior problem, the 
optimal gas density history in the exterior case becomes 
almost stationary only for d > 2. Furthermore, there is 
a subtle but important difference in the stationary MFT 


formulations of the interior and exterior problems, as we 
explain in the following. 

In the next Section we formulate the MFT in the con¬ 
text of the interior survival problem. In Sec. m we 
apply it to the non-interacting RWs in one dimension, 
where V is known, see Eqs. © and ©. In this case 
we can solve the full time-dependent problem exactly. 
The solution shows that the optimal density profile for 
this V is time-independent for most of the time. The 
full time-dependent solution will help us identify the cor¬ 
rect stationary formulation of the MFT problem. In Sec. 
m we apply this stationarity ansatz to an arbitrary in¬ 
teracting diffusive gas in any dimension. This yields a 
stationary equation which can be simplified further upon 
a transformation of variable. In Sec. 0 we apply this 
procedure to the SSEP. In subsection IV Al we present ex¬ 
act results for the SSEP in d = 1: for the stationary 
optimal density profile and for the survival probability. 
We verify these results, in the same Section, by solv¬ 
ing the full time-dependent MFT problem numerically. 
In subsection IV Bl we solve the stationary MFT problem 
for a rectangular domain, d = 2. In subsection IV Cl we 
study, analytically and numerically, the SSEP survival in 
spherical domains. In Sec. IV Dl we identify a universal 
behavior of the solution in the high-density limit of the 
SSEP: inside a domain of any shape and in any dimen¬ 
sion. Finally, in Sec. IVII we present a general solution of 
the gas survival problem in d = 1 which holds for a whole 
family of interacting diffusive gases. Section IVlTl presents 
a brief discussion of our results. For completeness, in Ap¬ 
pendix 0 we calculate the survival probability of a gas 
of RWs from the microscopic perspective. Appendix [B] 
extends the one-dimensional solution of the MFT equa¬ 
tions for the RWs, presented in Sec. cm to an arbitrary 
dimension and arbitrary domain shape. 


II. MACROSCOPIC FLUCTUATION THEORY 
OF PARTICLE SURVIVAL 


The starting point of the macroscopic fluctuation 
theory (MFT) [101 ] is fluctuational hydrodynamics: a 
Langevin equation for the fluctuating gas density g(x, t ): 


d t q 


V • [D(q)\/q] + V • 


Vo{q)ri(x,t) , 


( 8 ) 


where rj(x, t) is a zero-average Gaussian noise, delta- 
correlated both in space and in time. Equation © pro¬ 
vides an asymptotically correct large-scale description of 
fluctuations in a broad family of diffusive lattice gases 
0. At the level of fluctuational hydrodynamics, a dif¬ 
fusive gas is fully characterized by the diffusivity D(q) 
and additional coefficient, cr(q ), which comes from the 
shot noise and is equal to twice the gas mobility coeffi¬ 
cient [jjf ■ For example, for the non-interacting RWS one 
has D(q) = D 0 = const and a{q) = 2D 0 q , whereas for 
the SSEP D(q) = Dq = const and a(q) = 2D 0 q(l — q) 

EMI. 
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The MFT equations are essentially the saddle-point 
equations of the path-integral formulation, correspond¬ 
ing to the weak-noise limit of Eq. (0) (TO, 113, E5J . The 
MFT theory employs the typical number of particles in 
the relevant region of space as a large parameter. It al¬ 
lows to calculate the optimal path of the system: the 
most probable density history conditional on a specified 
large deviation. If the large deviation is described in 
terms of a spatial integral constraint, this constraint can 
be accommodated via the Lagrange multiplier formalism 
and provides a problem-specific boundary condition in 
time H3- 

Suppose we are interested in the probability that N 
particles are absorbed by the domain boundary by time 
T. (We will ultimately consider the limit of N = 0.) This 
defines an integral constraint on the solution: 


[ dx[n 0 - g(x,T)] = N. 
Jn 


(9) 


The same type of constraint appears in the exterior prob¬ 
lem (@, H] • A similar constraint also appears in the prob¬ 
lem of statistics of integrated current through a lattice 
site during a specified time [2i| - l28l . Referring the reader 
to Ref. [24] for a detailed derivation, we will only present 
here the resulting MFT equations and boundary condi¬ 
tions. 

The MFT equations can be written as two coupled 
partial differential equations for the optimal density field 
g(x,f) (the “coordinate”) and the conjugate “momen¬ 
tum” density field p(x, t): 

dtq = v ■ [D(q)Vq - a(q)S7p] , (10) 

d t p = -D(q)V 2 p-- cr'(q)(Vp) 2 , (11) 

where the prime denotes the derivative with respect to 
the single argument. Equations m and (EH) are Hamil¬ 
tonian, 


d t q = 6H/Sp, d t p = 

— SH/6q. 

(12) 

H[q(x, t),p(x, f)] = 

f dxT-L 

(13) 


Jn 


is the Hamiltonian, and 

n(q,p) = -D(q)S7q • Vp + \cr{q){\/pf (14) 

is the Hamiltonian density. The absorbing boundary im¬ 
poses zero boundary conditions for q: 

q(x € dCl,t) = 0. (15) 

Since the values of q are fixed at the boundary, the con¬ 
jugate field p must vanish there [1(3, l23l. [29j: 


(16) 


For the RWs and SSEP this boundarvcondition was de¬ 
rived from the microscopic models [23j, Ell- Although 
a general proof of Eq. © is unavailable [sol ], its va¬ 
lidity has been verified in many examples, see Refs. 

MEam. 

The boundary conditions in time are the following. For 
the density q we choose a deterministic initial condition 

q(x £ fl,t = 0) = n 0 . (17) 

The boundary condition in time for p(x, t = T) follows 
from the integral constraint (0), accounted for via a La¬ 
grange multiplier [24| : 

p(xe!l,( = T) = A, (18) 

p(x € cA2, t = T) = 0, 

where the Lagrange multiplier A is ultimately set by the 
constraint ©• The zero-absorption limit N = 0, that 
we are interested in here, corresponds to the limit of 
A -A +oo ! m. In this limit the particle flux to 
the boundary vanishes at all times 0 < t < T. 

Once the MFT equations with the proper boundary 
conditions are solved, we can calculate the action S that 
yields — In V up to a pre-exponential factor: 

— lnP ~ S' = f dt f dx(jpd t q — 'H) 

Jo Jn 

= \ [ dt [ dx <r(q) (Vp) 2 . (19) 

z Jo Jn 

III. MFT OF RANDOM WALKERS IN ONE 
DIMENSION: TIME-DEPENDENT SOLUTION 
AND STATIONARY ASYMPTOTIC 


A one-dimensional domain can be set to be an interval 
of length 2 R, centered at the origin. For the RWs the 
MFT equations become: 


d t q = D 0 d 2 q-2D 0 d x (qd x p), 

(20) 

d t p = —D 0 d x p — D 0 (d x p) 2 . 

(21) 

The boundary conditions in space are 

g(M =R,t) = 0, 

(22) 

p(|x| = R,t) = 0. 

(23) 

The boundary conditions in time are 

o 

£ 

II 

o' 

II 

''5< 

(24) 

p(x,t = T) = X6(R— |x|), 

(25) 


where 6 {...) is the Heaviside step function. 

As many other large deviation problems for the RWs, 
the problem (BUl) - E5l) is exactly soluble using the Hopf- 
Cole transformation Q = qe~ p and P = e p , defined by 
the generating functional 

/ R rR 

dx Fi ( q , Q)= dx [q In (q/Q) - q]. (26) 

-R J-R 


p(x € <9f2, t) = 0. 
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In the new variables the Hamiltonian density is 


H(q,p) = -D 0 d x Qd x P, (27) 


and the action can be expressed using only the initial and 
final states of the system (see the Appendix of Ref. 


-In V~S = 


dt / dxD 0 q(d x p) 


(28) 


rR 

/ dx[Q(PlnP-P+l)\ \ 7 . (29) 
J-R 


The transformed MFT equations are fully decoupled: 


d t Q = D 0 dlQ, (30) 

dtP = ~D 0 d 2 x P , (31) 

and the transformed boundary conditions are: 

Q{\x\ = R,t) = 0, (32) 

P(\ x \ = R,t) = 1, (33) 


and 


for the diffusion equation (1301) . The solution of the latter 
equation is 


Q{x,t) = n 0 



, G(x,x',t ) 

X l + (e A -l )g{x',T)’ 


(39) 


Transforming back to q and p, and taking the zero- 
absorption limit of A —> oo, we obtain: 

q(x,t) = n 0 g(x,T-t) ( dx' G ^ ’ 2 ’ ^ , (40) 

J-r 9{x',T) 

v(x,t) = d x p= d x kig(x,T -t), (41) 

r R 

— lnP ~ —no / dxlng(x,T). (42) 

J-R 


We are interested in the long-time limit, T R 2 /Dq. A 
close inspection of the exact relations (1551) and (l38l) re¬ 
veals an important feature that plays a crucial role in our 
further analysis. For D 0 t/R 2 1 and D 0 (T—t)/R 2 1, 

that is, outside of narrow boundary layers (in time) of 
typical width R 2 /Dq around t = 0 and t = T, the func¬ 
tions G(x, x' , t) and g(x, T — t) are well approximated by 
their lowest modes, n = 1 and n = 0, respectively: 


= = PiiJVoj- ‘ 34) 

P{x,t = T) = l + {e x -l)0(R-\x\). (35) 


The solution of the antidiffusion equation (l3lT) backward 
in time is obtained by integrating the “final” condition 
(1551) with Green’s function G(x,x',T — f), where 


G(x,x',t ) = 


sin 


7 vn(x + R) 


x sin 


n =1 

7T n(x' + R) 
2 R 


2 R 


^ 2 r» 2 P, 

g 4 R 2 


(36) 


The integration yields 

rR 

P(x,t) = 1 + (e A — 1) / dx'G(x,x',T — t) 

J-R 

= 1 + (e A — l)g(x, T — t), (37) 


where 


fR 

g{x,t) = / dx' G{x,x' ,t) 
J-R 


=E 

n=0 


7r(2n + 1) 


sm 


7r(2n + l)(x + R) 
2R 


2 (2n + l ) 2 P a t 

Tr 2 . 


(38) 


Evaluating P(x, t = 0) from Eq. (1571) and using Eq. (1551) . 
we obtain the initial condition, 


Q(x,t = 0) 


n 0 

1 + (e A — l)g(x, T) ’ 


G(x, x', t) 



(43) 


. m , 4 /7TXN K 2 D 0 (T-t) 

g(x, T -t) ~-cos[—J e ^ . (44) 

Plugging these approximations into Eqs. (150)) and SB, 
we obtain time-independent expressions: 


q{x) 

= 2no cos 2 ^ 

' 7TX \ 

2r) ’ 

(45) 

v{x) 

7T 


(46) 

=-tan 

2 R 

(w) 


That the density profile q(x , t ) is stationary most of the 
time is clearly seen in Figure [T| which shows the time- 
dependent solution at different times. Notice also that it 
is the momentum density gradient v(x , t ) = d x p , and not 
the momentum density p(x, t) itself, that stays almost 
stationary. This is in contrast to the exterior survival 
problem in d > 2. There the density profile is also almost 
stationary, but it is the momentum density, and not only 
its gradient, that is almost stationary Q. 

Using Eqs. da and (l42l) . we obtain the leading-order 
term of the long-time survival probability: 


— In V ~ no 


tt 2 A)T 
2 R 


(47) 


which is simply the action evaluated for the stationary 
solutions (1551) and (1551) on the entire interval 0 < t <T. 
This result agrees with Eqs. 0 and © as to be expected. 
Note that the stationary approximation remains accurate 
even when Dot/R 2 and Dq(T — t)/R 2 are of order unity. 
This is because the sub-leading terms in Eqs. m- m 
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x 



t 


where A is a yet unknown integration constant. Note 
that, for the exterior problem, a similar constant vanishes 
0- This mathematical difference between the internal 
and external problems is crucial. Although its effect on 
the solution is immediate, its physical interpretation is 
somewhat elusive. 

After plugging Eq. C51) into Eq. (15U1) we obtain a single 
second-order equation 

(VI 2 

q" - + 2Aq = 0 (51) 


which needs to be solved subject to the boundary condi¬ 
tion (1221) and a normalization condition, following from 
mass conservation: 


f R dxq(x) 


2 R 


= n 0 . 


(52) 


To solve Eq. m we make a transformation of variable 


q(x) = u 2 [x) 


(53) 


and obtain harmonic oscillator equation for u(x): 

u" + Au = 0. (54) 


FIG. 1: (Color online) The exact time-dependent optimal 
density profile (SSI) showing the formation of a stationary so¬ 
lution in do = 1, T = 10 and no = 0.5. Upper panel: the 
density profiles at times t = 0 (solid line), t = 0.001 (dashed 
line), t = 0.05 (dotted line), t = 1.5, 6 and 8.5 (three indis¬ 
tinguishable solid lines), and t = T = 10 (dash-dotted line). 
The lower panel shows q(x = 0, t) versus time. The narrow 
boundary layers at t = 0 and t = T are clearly seen. 


include a large factor 2n 2 in the exponent. This explains 
the narrowness of the boundary layers on Fig. [l] where 
D 0 T/R 2 = 10. 

Therefore, at sufficiently long times, the leading-order 
contribution to the survival probability of the RWs in¬ 
side a domain with an absorbing boundary comes from 
stationary solutions for q and v = d x p. The stationary 
solution for v solves the equation 

d t v = d x (- D 0 d x v - D 0 v 2 ) , (48) 

which is obtained by differentiating Eq. m in x. 

The stationarity of the optimal density profile, which 
gives the leading-order contribution to the survival prob¬ 
ability, is not unique to the one-dimensional case. We 
prove in Appendix [B] that this property holds for the 
RWs in any dimension and in an arbitrary domains. 

Now let us revisit the problem for RWs in one di¬ 
mension by directly looking for stationary solutions of 
Eqs. (l20l) and (|48jl . Setting the mass flux in Eq. (l20l) to 
zero, we end up with the following set of ordinary dif¬ 
ferential equations, which we will call stationary MFT 
equations: 


q{x) = 2 qv, (49) 

v'(x) = — v 2 — A, (50) 


The solution is u = B sin 


\/A (x + x 0 ) 


, where B , and 


xq are integration constants. The boundary conditions, 


zt(|x| = R,t) = 0, (55) 

set xq = R and A = m 2 ir 2 / (4 R 2 ), where m = 1, 2,.... 
Imposing the mass conservation (1521) . we obtain a family 
of solutions, parameterized by to: 


q m {x) = 2n 0 sin 2 


TO7r(x + R) 
2 R 


(56) 


The corresponding x-solutions, calculated from Eq fiiJl) . 
are the following: 


TO7T 

Vm\X) = —— cot 


2 R 


TO7r(x + R) 


2 R 


Now we can calculate the action from Eq. m- 


r R 


- In Vm ^ Sm = 


= m no- 


dt / dxD 0 q m v 2 
J-R 

it 2 _DoT 

2 R ' 


(57) 


(58) 


As we know from the full time-dependent solution, only 
the “fundamental mode”, to = 1 is selected by the ac¬ 
tual dynamics, see Eqs. C£D and m- Not surprisingly, 
this solution has the minimum action, see Eq. (1551) . We 
argue that the same feature (selection of the lowest sta¬ 
tionary mode) holds for the RWs in all spatial dimen¬ 
sions, see Appendix [Bj Furthermore, it also holds for a 
whole class of interacting diffusive gases. In the next 
section we derive stationary MFT equations for an arbi¬ 
trary diffusive gas. We then solve them for the SSEP, and 
support our findings by a numerical solution of the full 
time-dependent MFT equations for the SSEP in d = 1. 
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IV. STATIONARY MFT EQUATIONS FOR AN 
ARBITRARY DIFFUSIVE GAS IN ARBITRARY 
DIMENSION 


We start with the general MFT equations m and 
ED- Taking the gradient of Eq. ED. we obtain 


dtv = V 


-D{q)S7 v 



(59) 


where v = 'S/p. Now we look for time-independent solu¬ 
tions, g(x) and v(x) of Eqs. (fTUl) and (1551) . Equation (1551) 
yields 


V • v 


-±a'(g)v 2 - A 
D(q) 


(60) 


where A is an integration constant to be determined later. 
In its turn, Eq. ED yields a zero divergence of the mass 
flux, so that the mass flux is a solenoidal vector field. 
In the survival problem, this vector field must have zero 
normal component to the domain boundary. Using these 
two properties one can show (see Appendix A of Ref. 
mi) , that the minimum of the action is achieved when 
this vector field vanishes identically. Therefore, we arrive 
at the equation 


S7q 


^(g)v 

D(q) ■ 


(61) 


Essentially, this equation states that, for the optimal pro¬ 
file, the fluctuation contribution to the flux exactly coun¬ 
terbalances the deterministic flux. Expressing v from 
here and plugging it into Eq. ED. we obtain a closed 
equation for q: 


V • 


\Mq) 


S7q 


Ay/cr(g) 

D(q) 


(62) 


so that D(q) and cr(q) only enter through the combina¬ 
tion D/y/a. Introduce a transformation q = f(u ) that 
satisfies the equation 


D[f(u)] 
\/o [/(«)] 


/'(«) = I- 


(63) 


Performing the integration (and assuming that the inte¬ 
gral converges), we can define the function f{u) implic¬ 
itly: 


An additional constraint on u(x) comes from the mass 
conservation: 


In dx f ( u ) 

V 


n 0 , 


(67) 


where V is the domain’s volume. 

Nonlinear equations similar to Eq. (1651) appear in a 
host of physical problems. Probably the best known is 
the equation for stream function of an ideal incompress¬ 
ible fluid in two dimensions, where —Af(u) is the vortic- 
ity f32l - f3ll] . Although the original time-dependent MFT 
problem has a unique solution, the stationary problem, 
defined by Eqs. ED-ED. may have multiple solutions, 
and we will need to address the ensuing selection prob¬ 
lem. We already encountered this feature in the previous 
section, dealing with the RWs. 

The stationary formulation makes it possible to deduce 
the scaling of the optimal density profile, and of the sur¬ 
vival probability, with the linear system size L. We note 
that Eq. (1651) remains invariant upon rescaling of x by L , 
and of A by y/L. In the rescaled coordinates Eqs. (I66P 
and ED become 


u(x £ <9Q) = 0. 

fn dx f ( u ) . _ 


( 68 ) 

(69) 


where fl is the rescaled domain, and V is its volume. 
Extracting v from Eq. ED. substituting it in Eq. ED. 
and using Eq. (1631) , we obtain a simple expression for the 
action in terms of the new variable u: 

q 1 T d—2 r 

— ln'P~S' =- / dx[Vu(x)] 2 , (70) 

2 J h 

where u is the solution of the rescaled problem. As we 
can see, the L d ~ 2 scaling, previously observed for the 
RWs, see Eq. (ED, holds for a whole class of interacting 
gases. In particular, In V is independent of L for d = 2. 

We can also see how the survival probability depends 
on the diffusivity and mobility of the gas. Suppose that 
we can express D{q) and a(q) as D(q) = D 0 D(q) and 
a(q) = Do&{q), where Dq = D(no), and D , d are dimen¬ 
sionless functions of the dimensionless density q. Then, 
from Eq. (I6dl) we have u = y/Do F{q ), where F is a di¬ 
mensionless function of q determined solely by D and d. 
Using this relation in Eq. ED. we obtain 

-InV ~ S = D 0 TL d - 2 s{n 0 ), (71) 


[ f , D{z) 

/ dz r — = u. (64) 

Jo V a ( z ) 

As a result of this transformation, Eq. (l62ll becomes 

S/ 2 u + A f(u) = 0. (65) 

The boundary condition for u is the same as for q = f(u): 

u(x £ d Q) = 0. (66) 


where s(no) is a dimensionless function determined by 
the domain shape, and specialized for each model only 
via D and d. Comparing Eq. ED with Eq. 0 , we see 
that the particle interaction manifests itself only in the 
rescaled action s(n o). The same feature has been ob¬ 
served for the exterior survival problem [?]. 

For the RWs Eq. (1M1) yields / = u 2 /(2D 0 ), while 
Eq. (l65l) becomes the Helmholtz equation 

S7 2 u + (A /Do) u = 0, 


(72) 
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which admits analytical solutions for domains of simple 
shapes. In one dimension this equation coincides, up to 
a redefinition of A, with Eq. (l54l) of the previous section. 

Let us return to Eq. (1551) . For interacting diffusive 
gases the function f(u) is nonlinear. Still, Eq. (1551) can 
be solved analytically in one dimension, and we will ex¬ 
ploit this fact in Sec. eh In higher dimensions such a 
general solution is unavailable. Quite a few particular 
solutions, in different geometries, have been found for 
special choices of nonlinear /(it) (UHH]. Among them 
there is the case of f(u) ~ sin it, when Eq. (l65l) becomes a 
stationary sine-Gordon equation. Fortunately, this par¬ 
ticular case describes the well-known simple symmetric 
exclusion process (SSEP). As many other lattice gases, 
the SSEP behaves in its dilute limit as RWs, so that 
/(it —> 0) ~ it. Therefore, we will demand that the non¬ 
linear solution for the SSEP cross over at low densities to 
the (fundamental mode) of the Helmholtz equation (l72l) . 



FIG. 2: A plot of no(t') from Eq. (1821) . 


the coordinate rescaling, presented at the end of section 
hyi we can set L = R = 1, whereas Eqs. (1551) and m 
become 


V. SSEP: A STATIONARY SINE-GORDON 
EQUATION 

Substituting D{q) = 2D 0 q and cr{q) = 2D 0 q{l — q) 
into Eqs. (1541) and (l65l) . we arrive at the stationary sine- 
Gordon equation 


U{\x\ = 1) = 0, (78) 

/ dx sin 2 — 2no■ (79) 

(80) 


The general solution of Eq. m can be written as 


V 2 C/ + Csin U = 0, (73) 

where U = \]2/D 0 u , C = A /Do, and 

9 = f(u) = sin 2 (^j . (74) 

The survival probability is given by 

rpTcL— 2 r 

— In V~S = —-— / dx[Vu(x)] 2 
2 J n 

= D 0 TL d - 2 s(n 0 ). (75) 


where 


VC{x- Fxo/i'J j , ( 81 ) 

where sn(. .A is the Jacobi elliptic sine function, see 
e.g. Ref. [381 ]. whereas v and xo are integration con¬ 
stants. The boundary condition (1751) sets xq = 1 and 
C = 7ti 2 K 2 (^), where K(^) is the complete elliptic inte¬ 
gral of the first kind, and m = 1,2 .... The parameter v 
is uniquely determined by Eq. (EH) which gives, for any 
m, 


U(x ) =2 arcsin < \fv sn 


e <y) 

K(^) 


no, 


( 82 ) 


s(n 0 ) = 7 /rfx[Vtl(x)] 2 . (76) 

4 J n 

We will now solve Eq. (l73l) in several geometries. 


A. SSEP survival on an interval 

As we show in section IVI1 the one-dimensional case is 
exactly soluble in quadratures for any diffusive gas for 
which a stationary solution exists. In this section we find 
the explicit solution for the SSEP. For d = 1 Eq. (l73l) 
coincides with the equation of mathematical pendulum: 

cpU 

^L + CsmU = 0. (77) 

As for the RWs in section EH we set our domain to be 
a segment of length 2 R centered about the origin. With 


where E(^) is the complete elliptic integral of the second 
kind. The plot of no versus v is shown in Fig. [2] As one 
can see, there is a one-to-one mapping between 0 < no < 
1 and 0 < v < 1. 

What is left is to select the correct stationary so¬ 
lution out of the family of solutions parameterized by 
to = 1,2,.... We note that, in the dilute limit of the 
SSEP the solution must coincide with that for the RWs. 
This argument, and the action minimization, select the 
fundamental mode to = 1. Substituting U in (1741) . we 
arrive at the stationary ^-profile: 


q{x) 


VSV? [K(^)(a: + 1), v) = v 


cn 2 [K(u)x,u] 
dn 2 \K{v)x, v\ ’ 


where cn(...) and dn(...) are Jacobi elliptic functions. 
This solution is shown by the solid line in the upper panel 
of Fig. [3] 










0.8 
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FIG. 3: (Color online) Formation of a stationary optimal 
density profile of the SSEP in one dimension, see also Fig. 
[4] Upper panel: Solid line: stationary density profile (1831) . 
Symbols: density profiles found numerically for Do = 1 and 
t — 6, 12 and 18 (indistinguishable). Dotted line: numeri¬ 
cal profile at t = T = 24. Dashed line: the initial condition 
q(x,0) = no = 0.5. Lower panel: q(x = 0,1) found numeri¬ 
cally; the narrow boundary layers at t = 0 and t = T = 24 
are clearly seen. The analytical result for s from Eq. (1861) is 
s — 3.5137.... The numerical result from Eq. (1191) is 3.545. 


We can also calculate v(x) = d x p from the one¬ 
dimensional version of Eq. (ED). Going back to the phys¬ 
ical (non-rescaled) coordinate x, we obtain 


}{x) = -- 


K(u) sn 

K(GD) .. 

R 

dn 

i<M(d .. 

R ’ w 

Ren 

K(G(D .. 

R > w 



(84) 


Notice that v(x) cs =F (R — M) -1 as x —>■ ±R , in the same 
way as in the exterior survival problem [?]. The v(x)- 
profile from Eq. (1841) is shown, by solid line, in Fig. [4j 
Having found the stationary profile of q (or U), we can 
evaluate the survival probability from Eqs. CD and G0: 


— InD ~ 


D 0 T 
4 R . 
DpT 

R j ~ i 
D 0 Ts(u) 


dx (U x ) 2 


-l 

r 1 


dxK 2 (v) ncn 2 [K(n)(a: + 1), v\ 


(85) 



FIG. 4: (Color online) Formation of a stationary optimal pro¬ 
file of v(x, t) for the SSEP in one dimension, see also Fig. 
[3] Solid line: the stationary profile (1841) . Symbols: numeri¬ 
cal profiles at t = 4, 12 and 18 (indistinguishable). Dotted 
line: numerical profile at t=0. The parameters are Do = 1, 
no = 0.5, and T = 24. 


where 


s(is) 


2K 2 (u) 


' EM 

KM 


+ u-l 


( 86 ) 


Equations (1521) and (1551) determine the rescaled action 
s = s(no) hr a parametric form. The low- and high- 
density asymptotics of s(no) are the following: 


s(n 0 ) 


2 ’ 

2 

1 - no ’ 


no < 1, 

1 — n 0 <C 1, 


(87) 

( 88 ) 


The low-density asymptotic coincides with that for the 
RWs, see Eqs. © and ©. The high-density asymptotic 
formally diverges as no approaches 1. This divergence, 
however, is cured when no becomes very close to 1, as 
explained below in this subsection. Figure [5] shows a plot 
of s(n 0 ), alongside with the asymptotics (1571) and (1551) . 

Also shown in Fig. [Dare numerical results obtained by 
solving the full time-dependent MFT equations m and 
CD for the SSEP in one dimension, using the boundary 
conditions in time CD and m- The numerical solution 
was obtained with a modified version of the iteration al¬ 
gorithm used in Ref. [jij for the exterior problem. Figures 
[3] and 0 show the time-dependent numerical solutions for 
q and v respectively, at different times. Apart from nar¬ 
row boundary layers at t = 0 and t = T, the solutions are 
very close to the analytical stationary solution. The nu¬ 
merically evaluated rescaled action s(n o), shown in Fig. 
0 is also in very good agreement with the analytical re¬ 
sults. 

Now let us return to the high-density limit where the 
MFT action (1551) tends to diverge. As 1 — no -C 1, or 
1 — v <C 1, we can approximate snz = tanh z + 0(1 — v), 
and K(j/) ~ (1 — no) -1 . As a result, 


U(x) ~ 2 arcsin 


tanh 


1 - M 

5 


R 


(89) 
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FIG. 5: (Color online). Solid line: the rescaled action s(no) 
from Eqs. (1821) and (1861) . Dashed line: the RW asymptotic 
(|57|). Dotted line: the high-density asymptotic (1881) . Sym¬ 
bols: s(no) computed from the numerical solution of the time- 
dependent MFT equations. 

where S = 1 — no <C 1. The resulting density profile 
q = sin 2 ( U/ 2) describes two kinks, with characteristic 
width 5 = 1 — no <C 1, located close to the ends of the 
interval: 


in domains of different shapes. An apparent divergence 
of S at no —> 1 also appears there (see below), and it is 
also cured at the microscopic level. 


B. SSEP survival inside a rectangle 

Here we will solve the stationary sine-Gordon equation 
(1731) inside a rectangular domain with dimensions L x and 
L y . After rescaling the coordinates, the dimensions of 
the rectangle become 1 and A = L y /L x , in the x and 
y directions, respectively. The equation must be solved 
with zero boundary conditions, whereas Eq. (El reads 



(91) 


As we explain shortly, an infinite family of solutions to 
this problem can be obtained by the method of “general¬ 
ized separation of variables” [34J. In the dilute limit, one 
of this solution coincides with the fundamental mode of 
the Helmholtz’s equation (El, 

U = 4,/no sin ( nx ) sin , (92) 

obeying the zero boundary conditions and the normaliza¬ 
tion condition m, where one should replace sin([//2) by 
U/2. We argue, therefore, that this solution yields the 
true stationary optimal density profile. 

The generalized separation of variables employs the 
ansatz 


q(x) ~ tanh 2 • 


(90) 


The action mostly comes from the kinks, and each kink 
contributes — 1/S to the action, leading to the asymp¬ 
totic (1551) . Now we will see how the apparent diver¬ 
gence of the action (1551) at % -f 1 is cured. The 
MFT is only expected to apply when the length scales 
that it describes are much greater than the lattice con¬ 
stant a. Restoring all units, we can express the kink 
width as S x R = (1 — ano)R. The MFT is valid 
when this quantity is much greater than a, that is when 
1 — ano ^ a/R. On the other hand, exactly at close 
packing, no = 1/a, the survival probability of the SSEP 
is equal to the product of probabilities Pi .2 that each of 
the particles adjacent to the boundary does not hit the 
boundary during time T. Each of these probabilities is 
Pi ,2 = exp(— D 0 T/a 2 ), so the total survival probability is 
equal to exp(— 2D 0 T/a 2 ) and is of course finite. As one 
can see, the crossover between the macroscopic result, 
InP mft — —2DqT/{ ciR(l — ano)], and the microscopic 
result, InP = —2 D 0 T/a 2 occurs at 1 — ano ~ a/R, when 
the MFT theory breaks down. 

As we will see in the following sections, the kink solu¬ 
tion (15U1) plays an important role in the high-density be¬ 
havior of the stationary solutions in higher dimensions, 


U(x, y) = 4 arctan [f(x)g(y)] (93) 

which yields two uncoupled equations for / and g (see 
Ref. [34| for a detailed derivation): 

iff = nf 4 + mf 2 + k, (94) 

(. g'f = kg 4 - (m + C)g 2 + n, (95) 

where to, n, and k are arbitrary parameters. Each of 
Eqs. (l94]l and (l95ll describe conservation of energy of an 
effective classical particle in a potential. As one can see, 
the particle motion is confined, and the resulting solution 
exhibits the correct low-density asymptotic (1551) . if and 
only if — C < m < 0, n > 0, k > 0, and (to + C) 2 > 4 kn. 
In this regime of parameters the solutions for / and g are 
elliptic functions [38j : 


/ 

g 


I —mv\ 
n( vi + 1) 

I (to + C)v 2 
k(v 2 + 1) 


sn 


^i + l 


(x + ci),vi 


sn 


to + C 


V2 


1 


(y + c2),v 2 


(96) 

(97) 


where c\ and C 2 are the integration constants of the first 
order equations (1M1) and (1^51) . and the constants v\ and 
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1*2 are given by m, n, and k via the relations 


Vi 

kn 

(98) 

(l + ^l) 2 " 

9 ’ 

m z 

^2 

kn 

(99) 

(1 + ^2) 2 

(to + C ) 2 ’ 

Imposing the zero boundary 

condition we obtain 


/ —m 

V (1 + ZT) ~ 

2?7ZiK(z/ 1 ), 

(100) 

m + C 

2m 2 K(p 2 ) 

(101) 

V (1 + Vi) 

A 

Cl =C 2 = 

0, 

(102) 


where mi and m 2 are positive integers. Similarly to the 
one-dimensional case, we must put mi = m 2 = 1. Now 
we can solve Eqs. (1951) - (I101I) for zq and obtain an expres¬ 
sion for fg in terms of zq alone: 


fg = (zqz* 2 ) 1/4 sn [2K(zq)x, zq] sn 


2K (u 2 )y 




, (103) 


where zq and zq are related by the equation 

K(zq) 4 zq = [K(zq)A] 4 zq. (104) 

Now we use mass conservation where we substitute 


q = sin 2 



4 (fg) 2 

[1 + (/ 5 ) 2 ] 2 ' 


(105) 


Thus all the constants are determined implicitly. The 
survival probability is given by Eq. ED, which we can 
rewrite as 


-In V~S = 


DnT 


dx 


dy (d x U) 2 + (d v U) 2 


= D 0 Ts(n 0 , A). 


(106) 


As to be expected from Eq. ED- the resulting probability 
is independent of the system size, but it strongly depends 
on the gas density. Figure [6] shows a two-dimensional plot 
of q(x,y) for no = 0.67 and A = 0.7. Figure 0 depicts 
s(no, A = 1) versus no, alongside with the low- and high- 
density asymptotics that we will now discuss. 


1. Dilute gas, no < 1 

In the dilute limit, no 1, our results coincide with 
those for the RWs. Indeed, by virtue of Eq. ED , u 
must be much smaller than 1 in this limit. Therefore, fg 
must be also small, and we can approximate U ~ 4 fg, 
and q cs 4 f 2 g 2 . Now, from Eq. ED, there must be 
zq, zq < 1. Therefore, as in the one-dimensional case, 
we can put sn(...,z*) ~ sin(...) and K(z*) ~ ir/2 in 
the expression for U. Further, Eq. (11041) yields in this 



x 


FIG. 6: (Color online) The stationary optimal density profile 
q(x,y) from Eq. (11051) . corresponding to the SSEP survival in 
a rectangle with aspect ratio A = 0.7 and no — 0.67. 



FIG. 7: (Color online). Solid line: the rescaled action for a 
square vs. no from Eqs. ED and (11061) with A = 1. Dashed 
line: the RW asymptotic (11081) . Dotted line: the high-density 
asymptotic JHED- 


limit z *2 — zqA 4 . The remaining constant zq ~ no/ A 2 is 
found from Eq. ED- As a result, we arrive at the cor¬ 
rect asymptotic (l92l) . The optimal density profile in the 
dilute limit, back in the original coordinates, is 


q ~ 4rio sin 


• 2 (irx\ . 2 (ny 


\L X 


sm 


(107) 


The survival probability is given by Eq. (11061) with 

s(n 0 , A) = 7r 2 ?r 0 ( A + . (108) 


The survival probability is maximum when the rectangle 
is a square. 
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2. Near close packing, 1 - r*o < 1 


C. SSEP survival inside a sphere 


When no is close to 1, the gas density is close to 
1 everywhere except in narrow boundary layers of the 
size 6 = 1 — no along the boundary. From Eq. (I105D . 
fg is close to 1. As v\ and v- 2 are also close to 1, 
K(zq) and K/tq) diverge as no 1. In this limit 
Eq. (11041) yields K{y 2 ) — K(iq)A. The explicit no¬ 
dependence can be obtained with the help of Eq. (EH): 
2K(m) ~ (1 + A)(A<5)-\ and 2K(n 2 ) ~ (1 + A)(<5) _1 . 
Using the asymptotic snz ~ tanh z, we obtain in the 
leading order: 


Here the stationary optimal density profile depends 
only on the radial coordinate r, and Eq. m becomes 

V 2 ,U + CsmU = 0, (113) 

where V 2 [/(r) = U rr + (d — l)U r /r is the radial Lapla- 
cian in d dimensions. Upon the coordinate rescaling 
r —> r/R , we need to solve the stationary sine-Gordon 
equation (I113D inside a sphere of unit radius. The bound¬ 
ary conditions, and the normalization condition, are 


fg 


tanh 


(1 + A)a; 


AS 


tanh 


(1 +A)y 


A<5 


tanh 

'(1+A)U 

AS 

+ tanh 

Rl+Ah/I 

A<5 

- 1 

tanh 

(l+A )(x+y) 
AS 



(109) 


This asymptotic is valid for x <1/2 and y < A/2; it can 
be extended to the other three quarters of the rectangle 
by obvious reflections. Away from the domain corners we 
can replace tanh(...) in the denominator of Eq. (11091) by 
unity. Using the resulting “approximate product rule” 
in Eq. (11051) . we obtain, after some algebra, the “kink” 
asymptotic of q(x, y) away from the domain corners: 


q — tanh 2 


'2(1 + A)x 

tanh 2 

I2(l + A)y] 

AS 


AS 


( 110 ) 


where x < 1/2 and y < A/2, and reflected formulas 
in the other three rectangle quarters. These asymp¬ 
totics describe kinks with the characteristic width £ = 
A<5/[2(1 + A)]. By analogy with one dimension, the ac¬ 
tion per unit length along the boundary is, in the leading 
order, \/£. Multiplying this expression by the perimeter 
2(1 + A), we extract the asymptotics 


U'(r = 0) = U(r = 1) =0, (114) 



(115) 

Then, from Eqs. (EH and CZHD, we obtain — In V ~ S = 
D 0 TR d ~ 2 s(no), where 


s(n 0 ) = ~^J drrd 1 ( u r) 2 , (H6) 

and fid is the surface area of the d-dimensional unit 
sphere. In the absence of general analytic solution of 
Eq. (11131) for d > 1, one can solve this equation numer¬ 
ically, and also explore analytically the low- and high- 
density limits. 

The first-order term of the density expansion of s(no) 
corresponds to the RWs, see Appendix iBl The next, rig- 
term can be obtained by treating the q 2 term of the MFT 
Hamiltonian of the SSEP perturbatively, similarly to how 
it was done in the exterior problem |7j. For example, for 
d = 3 the resulting correction is 

<5s(n 0 ) = 4 tt / drr 2 4 w (r)t£ w (r), (117) 

Jo 


s ( n o, A) ~ i- no + A + 2 ) ' 

Again, the survival probability is maximum when the 
rectangle is a square. 


where <7 rwM and urw(>') are the stationary optimal pro¬ 
files for the RWs, given by Eqs. (IB24I) and (IB25I) of Ap¬ 
pendix B. Evaluating the integral, we obtain Ss(no) = 
a rig, where 

o 4 

a = —— [2 Si (27t) — Si (47r)] = 38.7945 ..., 


3. Very long rectangle, A 1 


and 


Here, sufficiently far from the edges y = 0 and y = A, 
U(x,y) is almost independent of y, and close to the one¬ 
dimensional solution (ISTll . Therefore, the rescaled action 
s(n-o, A) is approximately equal to 

s(n 0 , A > 1) ~ 2Asi d (n 0 ), (112) 

where sid(no) is described by Eqs. (1821) and (l86l) . The 
factor A is due to additional integration along y, and the 
factor 2 appears because the one-dimensional result (l86l) 
was obtained for a segment of length 2, not 1. 



is the sine integral function. The resulting low-density 
asymptotic, for d = 3, is 

47T 3 

s ( n o) - — no + anQ. (118) 

Near close packing, 1 — no <C 1, the gas density q(r) 
drops from a value close to 1 to zero in a narrow boundary 
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FIG. 8: (Color online). Stationary optimal density profiles 
for the SSEP survival inside a sphere, d = 3. Solid lines: 
profiles, obtained by numerically solving Eq. (fTTifl) with three 
different values of C, corresponding to three different values 
of no, and using Eq. (ITU) . Curve 1 : C = 10.4 and no ~ 0.027. 
Dotted line: the asymptotic profile (IB24I) . Curve 2: C = 20 
and no — 0.31. Curve 3: C = 300 and no — 0.82. Dashed 
line: the asymptotic profile CEIIl- 


layer, of width 0(5) near the sphere r = 1. Correspond¬ 
ingly, U(r) =2 arcsin \Jq(r) rapidly drops from a value 
close to 7T to zero. As a result, we can neglect the first- 
derivative term in the radial Laplacian. This brings us 
back to the equation U rr + C sin U = 0 that we consid¬ 
ered in Sec. IV Al The boundary conditions (11141) are also 
the same as in the one-dimensional case, see Eq. eg. 
The only difference is in the normalization condition, 
Eq. (11151) which introduces the factor d. As a result, 


U(r) cs 2 arcsin 



d( 1 — r) 
1 - no 


and 



d( 1 — r) 
1 - n 0 


(119) 


( 120 ) 


Using Eq. (11191) . we obtain the high-density asymptotic 


s(n 0 ) = ^ [ drr d 1 (U r ) 2 ~ ^ , (121) 

4 Jo 1 - n 0 

In particular, for d = 3, 

127T 

S 3 d(no) — - -, 1 — n 0 <Cl. (122) 

1 - n 0 

For an arbitrary density no Eq. (11131) can be solved 
numerically: either by the shooting method or by artifi¬ 
cial relaxation. Examples of such solutions for d = 3 are 
shown in Fig. [8] Figure [9] shows the numerically found 
S 3 d(no), alongside with the asymptotic (11181) . its linear 
term only, and asymptotic (11221) . 



FIG. 9: (Color online). Solid line: the rescaled action s(no) 
for a sphere, d = 3, vs. no, obtained by numerically solv¬ 
ing Eq. ( 11131) and using Eq. (11161) . Also shown are the low- 
density asymptotic (11181) (dash-dotted line) and its linear 
term (dashed line). Dotted line: the high-density asymp¬ 
totics (11221) . 


field q stays very close to 1 across most of the domain, 
and drops to 0 in a narrow boundary layer of charac¬ 
teristic width 6 = 1 — no along the domain boundary. 
The density derivatives in the directions parallel to the 
boundary are negligibly small compared to the density 
derivative across the boundary. Therefore, we can ap¬ 
proximate the Laplacian in Eq. by a one-dimensional 
one, d 2 U/dr\ : where rj_ is a local coordinate normal to 
the domain boundary. As a result, the problem becomes 
effectively one-dimensional, and the solution of Eq. GSD 

is a one-dimensional kink, g(x) ~ tanh 2 (^/C r ±^ (we 

set r _l to vanish at the boundary). The action, Eq. (1751) . 
mostly comes from the boundary layer. In the leading 
order, the action per unit surface across the boundary is 
equal to VC. The total action is therefore A\/C, where 
A is the normalized domain’s surface area. \[C is deter¬ 
mined by the mass conservation: \[C ~ A/(V8), where 
V is the normalized domain’s volume. The final result is 


-\nV ~ S ~ D 0 TL 


d—2 


A 2 DnT 


V (1 — no) V (1 — ng) 


(123) 


This expression is in agreement with our high-density 
results (1551) . (I111D . and (11211) . 


D. SSEP survival in arbitrary domains near close 
packing 

Consider a domain of arbitrary shape, in any dimen¬ 
sion. As n o approaches 1, the stationary optimal density 


VI. SURVIVAL OF A GENERAL DIFFUSIVE 
GAS ON AN INTERVAL 

For a general diffusive gas in one dimension the sta¬ 
tionary density profile and the survival probability can 
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be found in quadrature. Indeed, for arbitrary D(q) and 
a(q) in one dimension, Eq. (1551) reads 


u" + A f'(u) = 0, (124) 

where f(u) is defined by Eq. (l64l) . By virtue of the scaling 
properties of the problem, it suffices to solve Eq. (112411 on 
the interval |a;| < 1, with the same additional conditions 
for u as stated in Eqs. © and ([79]) for U . Equation (11241) 
describes the motion of an effective classical particle with 
unit mass (u is the “particle coordinate”, x is “time”) in 
the potential V(u) = A/(it). Let us denote the particle 
energy by Av. Energy conservation yields a first-order 
equation: 


u = ±y/2A [v- f{u)} (125) 


Solving it with the boundary condition u[x 
we obtain: 


dh 


Jo \Jv- f{h) 


= V2A(x + 1). 


- 1 ) = 0 , 
(126) 


where we have again used f\q(x)dx = 2no . Using 
Eq. (11281) . we can rewrite Eq. (11321) as 

- hrU ~s(n°), (133) 

where the rescaled action s = s(no) is given in a para¬ 
metric form by the equation 


s(u) = 


dqD(q) 


o \/{y ~ q) cr(q) 


\v - n 0 (is)] (134) 


and Eq. (11301) . When specialized to the SSEP, Eqs. (11271) 
and (11341) yields (1531) and (1551) respectively. 

With these general results at hand, we can investigate 
the survival properties of a whole class of diffusive gases 
with known D(q) and cr(q) in one dimension: on the con¬ 
dition that the integral in the denominator of Eq. (11301) 
converges. 


VII. CONCLUSIONS AND DISCUSSION 


Changing the integration variable to z = f{h) and using 
Eq. (1551) . we obtain 


V D(2) 


Jo y/(v - z)a(z) 


= V2A{x + l). 


(127) 


Now we demand that q = v at x = 0, and express A 
through u: 


A=i 

2 


dz D(z) 


-i 2 


(128) 


Jo \J(v - z)a(z) 

An additional condition comes from mass conservation: 

f 1 [ v dx 

2no = / q(x)dx = 2 q —— (q) dq 

J -1 Jo dq 

qD(q ) 


= 2 dq 


o \/2A (v-q)a(q) 


(129) 


Equations (I128D and (I129D yield the dependence of no on 
v. 


<?£>(<?) 


n 0 (i/) = 


r dq n -- 

fodq- 


D(q) 


(130) 


>0 

Now we can evaluate the survival probability: 


T 


( du 


— In V ~ —— / dx l — 


2 R 


\dx 


(131) 


By virtue of Eq. (11251) and the definition q = f(u ), we 
obtain 

, -T, A T f\ f 2A T(u-n 0 ) , 1on , 

J \y-q{ x )\dx= ---, (132) 


We dealt in this work with the survival of a gas of inter¬ 
acting diffusive particles inside a domain with absorbing 
boundary. Employing the MFT formalism, we evaluated 
the long-time survival probability of the gas and its op¬ 
timal density history conditional on the survival. We 
found that this optimal density history is stationary dur¬ 
ing most of the process. As a consequence, the survival 
probability decays exponentially in time: inside domains 
of any shape in all dimensions. As we showed, the solu¬ 
tion of the long-time survival problem reduces to solving 
a nonlinear Poisson equation, where the nonlinear term 
is determined by D(q) and cr(q). In one dimension, this 
problem is soluble exactly for a whole class of diffusive 
gases. For the SSEP the nonlinear Poisson equation be¬ 
comes a stationary sine-Gordon equation, and we solved 
it in different geometries and dimensions. 

The dilute limit of the SSEP corresponds to non¬ 
interacting random walkers (RWs), where the problem 
reduces to finding the lowest positive eigenvalue /xi of 
the Laplace’s operator inside the domain Q, see Eq. ©. 
Near close packing the problem becomes effectively one- 
dimensional and can be approximately solved for any do¬ 
main shape and any dimension. What is the optimal 
domain shape, for a fixed number of particles and fixed 
volume of the domain, that maximizes the chances of 
long-time survival? Interestingly, both in the dilute limit 
of the SSEP, and near close packing, the optimal domain 
shape is the sphere. Indeed, in the dilute limit the mini¬ 
mum value of /ii is achieved for the sphere, as guaranteed 
by the Rayleigh-Faber-Krahn theorem (3§, [4y|. In its 
turn, near the close packing, the sphere is the minimizer 
of the surface area A at fixed volume, see Eq. (I123D . A 
natural conjecture is that the sphere maximizes the sur¬ 
vival probability of the SSEP at any gas density, but we 
do not have a proof. For other diffusive gases we do not 
know the minimizing domain shape. 
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It would be interesting to apply our approach to the 
“narrow escape problem”, where there is a small hole 
in the reflecting boundary of the domain. The survival 
probability [Tl] and the mean escape time [42j of a single 
RW in this system have been extensively studied. The 
MFT formalism can give an interesting insight into fluc¬ 
tuations in the escape of interacting particles. 

Another interesting extension would address the full 
absorption statistics: evaluating the probability that a 
specified number of particles are absorbed by time T. 
An exterior variant of this problem has been recently 
considered, for the SSEP, in Ref. @|. 

Finally, we emphasize that stationarity of the optimal 
gas density profile is a major simplifying factor in the 
large-deviation problem we have considered here. Cases 
of time-independence of the optimal gas density history 
are also encountered in other large-deviation settings in 
lattice gases [ 3 , :j§, , and they are intimately related 
to the “additivity principle” put forward by Bodineau 
and Derrida [45] . 
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at x^. This probability distribution obeys the diffusion 
equation (2H3 : 


d tP i = D 0 W 2 p i (A3) 

with the absorbing boundary condition: 

pi(x £ dfl, t) = 0. (A4) 

The initial condition is 

pi(x, t = 0) = S d (x - Xj), (A5) 

where S d is the d-dimensional Dirac delta-function. The 
solution to Eq. (IA3I) is the Green’s function of the diffu¬ 
sion equation: 


/9i(x,t,Xi) = G(x,x,,t) = ^ 4 , n(x)T n (x,)e ^ Dot . 

n= 1 

(A8) 

Here and p n are the normalized eigenfunctions and 
eigenvalues of the eigenvalue problem V 2 u + p 2 u — 0 
for the Laplace’s operator inside the domain with the 
boundary condition u(x £ dfl) = 0. We order the eigen¬ 
values by their magnitude pi < P 2 <Us < ■ ■ ■ ■ Plugging 
Eq. (IA6I) into Eq. (IA2I) one obtains Q: 


g(xi,T) = / dxG(x,Xi,t) 

J n 

OO r. 

= V^(x,:) / dx4/„(x)e-^ D ° T (A7) 

■n-l J n 


Appendix A: Survival probability of Random 
Walkers from a microscopic perspective 


2. Evaluating the sum in Eq. (ED 


For completeness, we present here a brief microscopic 
derivation of the survival probability of a gas of non¬ 
interacting Random Walkers (RWs) inside a domain S2. 
This quantity can be obtained from a single-particle sur¬ 
vival probability: 

— In V (T) RW = — ^2 In <7( X j, T) , (Al) 

i 

where g(xi,T) is the survival probability up to time T of 
a particle initially positioned at Xj , and the sum is over 
all particles. Therefore, one needs to calculate g(x,,T) 
and perform the summation. 

1. Calculating (?(xi, T) 

The single-particle survival probability g(x.i,T) is 
given by the expression 

g(xi,T) = [ dxpi(x,T, Xj), (A2) 

Jn 

where pi(x, t,Xj) is the probability distribution of the 
particle position, given its (deterministic) initial position 


When the number of particles in the domain Q is very 
large, the sum in ED can be approximated by the inte¬ 
gral: 

- lnU(T) RW ~ -n 0 f dx' In [ ff (x, T)], (A8) 

Jn 

Furthermore, at times much longer than the character¬ 
istic diffusion time, the infinite sum in Eq. ED can be 
approximated by its first term: 

<?(x, T) ~ 4'i(x) [ dx / 'F 1 (x / )e _/i i D ° r . (A9) 

Jn 

Plugging this approximation into Eq. (|A8[) . we obtain 
the long-time asymptotic of the survival probability pre¬ 
sented in Eq. ©. 

Appendix B: Solving the MFT equations for the 
RWs in higher dimensions 

The one-dimensional solution, presented in section Hill 
can be generalized to any simply-connected domain in 
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arbitrary dimension. Consider the MFT equations m 
and (flTT) for the RWs: 

d t q = V • [Do'Vq - 2 D 0 q'\/p\ , (Bl) 

dtp = -D 0 V 2 p-D 0 (Vp) 2 , (B2) 

The Hamiltonian density is 

n(q,p) =-D 0 \7q-Vp + D oq (\7pf . (B3) 

The absorbing boundary conditions are described by 
Eq. (1151) and (1161) . The boundary conditions in time are 
given by Eqs. (fTTl) and (fl8l) . 

As in one dimension, we solve the problem, using the 
Hopf-Cole transformation Q = qe~ p and P = e p , with 
the generating functional 


The resulting solution of Eq. (IB8j) is 


Q(x, t) =n 0 dx' 
Jn 


G(x,x',t) 


1 + (e A — l)g(x',T)' 


(B15) 


Now we calculate the action using Eq. (jB6|) . After some 
algebra, and taking the zero-absorption limit of A —> oo, 
we arrive at Eq. (IA8I) , which describes the continuum ap¬ 
proximation of the exact microscopic result ED- Trans¬ 
forming back to q and p, and taking the limit of A —> oo, 
we obtain: 


g(x,t) = n 0 g(x,T-t) f dx 1 X ’ ^ , (B16) 
Jn 5 (x',T) 

v(x,i) = Vp = Vln 5 (x,T-t). (B17) 


[ dxF 1 (q,Q) = ( dx [q \n(q/Q) — q] . (B4) 

J Q J Q 


/f2 JO, 

In the new variables the Hamiltonian density is 

H(q,p) = -D 0 VQ-VP, 


(B5) 


and, again, the action can be expressed through the ini¬ 
tial and final states of the system: 

— lnT^Rw — S= f dt f dxD 0 q(yp) 2 
do Jn 

= [ dx[Q (PinP — P + 1)] 

Jn 


The transformed MFT equations are decoupled: 

d t Q = DqV 2 Q } 
d t P = -D 0 V 2 P. 


The transformed boundary conditions, in space and in 
time, are: 


Being interested in long times, we observe that, outside 
the boundary layers of width L 2 /Do around t = 0 and 
t = T, one can approximate expressions Eq. E3 and 
m by the first terms of the corresponding series: 

G(x, x', t) ~ 4'i(x)4'i(x / )e-^ Dot , (B18) 

g(x,T-t)~\ Et(x) [ dx'^ 1 (x')e-rt Do{T - t) .( B19) 
Jn 


(B6) 

This approximation yields the stationary solution 


(B7) 

q(x) = n 0 V^l(x) 

(B20) 


, , „ V^i(x) 

v(x) = V P = #i(x) . 

(B21) 

(B8) 



(B9) 

whereas — In V is given by Eq. © . 



When fl is a circle of radius R (d = 2), we obtain 


Q(x e dfl, t) = 0, 

(BIO) 

q{r) = 

n oJo ( R ) 

JH z i) 

(B22) 

P(x £ dfl, t ) = 1, 

O( x ,i = 0) = p(xt = 0 ). 

(Bll) 

(B12) 

v(r) = 

Rj^m ’ 

(B23) 


and 


P(x,t = T) = 


e . 

1. 


x s n 

x € dfl 


(B13) 


where Jo(z) and J\(z) are Bessel functions, and z\ = 
2.4048 ... is the first positive root of Jo(z). The survival 
probability is described by Eqs. © and ©• 

When H is a sphere of radius R (d = 3), the stationary 
solution is 


Solving the anti-diffusion equation (IB9I) . we obtain 

P(x,t) = l + (e x -l)g(x,T~t), (B14) 

where g is defined in Eq. (lA7l) . Evaluating P(x,t = 0), 
we obtain the initial condition for the diffusion equa¬ 
tion 


q(r) = 
v(r) = 


2 noR 2 sin 2 (^) 
3? ' 


7T 

/nr \ 

1 

— cot 

— 

— _ 

[R 

kR) 

r 


(B24) 

(B25) 


Q(x,t = 0) = 


n 0 


1 + (e A — l)<?(x, T)' 


and the survival probability is described by Eqs. Q and 

©. 
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